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Abstract 

Background: This paper provides some clarifications regarding the use of model-fitting methods of kinetic analysis 
for estimating the activation energy of a process, in response to some results recently published in Chemistry 
Central journal. 

Findings: The model fitting methods of Arrhenius and Savata are used to determine the activation energy of a 
single simulated curve. It is shown that most kinetic models correctly fit the data, each providing a different value 
for the activation energy. Therefore it is not really possible to determine the correct activation energy from a single 
non-isothermal curve. On the other hand, when a set of curves are recorded under different heating schedules are 
used, the correct kinetic parameters can be clearly discerned. 

Conclusions: Here, it is shown that the activation energy and the kinetic model cannot be unambiguously 
determined from a single experimental curve recorded under non isothermal conditions. Thus, the use of a set of 
curves recorded under different heating schedules is mandatory if model-fitting methods are employed. 



Findings 

In a paper recently published in the Chemistry Central 
Journal, a kinetic study of the thermal decomposition of 
both aged and non-aged commercial cellulosic paper 
was presented, and the apparent activation energy (Ea) 
of the degradation reaction was determined for each 
case [1]. According to the authors, the Ea of the process 
is related to the breakdown of cellulose chains and, since 
the apparent activation energy of the process is found to 
decrease with the aging time of the cellulose paper, it is 
proposed that such evolution could be used to construct 
archaeometric curves. Three different model-fitting 
methods were used to determine the activation energy: 
the differential Arrhenius method, the integral Savata 
method and the Wyden-Widmann method. Also, the 
authors speculate with the possibility that the kinetic 
method selected influences the obtained Ea. Actually, 
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when using the Wyden-Widmann method, it is observed 
that only a limited number of data points around the 
DTG peak should be employed or else the Ea 
obtained would not fit that obtained by the other kin- 
etic methods. Finally, it is concluded that a first order 
model is the most suitable for describing the cellulose 
decomposition reaction. However, recent works have 
reached to different conclusions, finding a chain scis- 
sion model to be far more appropriate [2]. Such dis- 
crepancy is due to some fundamental misconceptions 
in the manner the kinetic methods are employed in 
Marini's work. Firstly, the activation energy for every 
sample studied was obtained by means of applying a 
model-fitting method to experimental data proceeding 
from a single non-isothermal run. Secondly, only the 
fit to two kinetic models, Fl and A2, were explored 
in the analysis. Basically, model-fitting methods of 
kinetic analysis consist of fitting the experimental data 
to a series of theoretical kinetic models, which are 
algebraic functions that reflect the relationship between 
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reaction rate and degree of conversion and can be 
related to the reaction mechanism. The model provid- 
ing the best linear fit is usually regarded as the cor- 
rect one, and the activation energy is deduced from 
the slope of the fit. Unfortunately, it has been long 
established that the activation energy cannot be reli- 
ably determined from a single non isothermal curve 
because the experimental data almost always provides 
a reasonable fit regardless the kinetic model selected 
[3-5]. Despite that significant flaw, such inappropriate 
practice is still nevertheless widely used. As a result, 
it is common that «th order models are incorrectly 
selected because they are often tested as the first op- 
tion for simplicity and a good fit is usually obtained. 
Here, we attempt to throw some light on the use of 
model-fitting methods and clarify such still wide- 
spread misuses. 



Results and Discussion 

Figure 1 includes a simulated oc-T curve, constructed as- 
suming a heating rate of 10 K min" and the following 
kinetic parameters: E a =150 kj/mol, A=10 s~ and a Fl 
(first order) kinetic model. The model-fitting methods of 
Arrhenius and Savata, those used in Marini's work 
[1], were selected to determine the activation energy 
of the simulated curve. Thus, from the simulated 
data, the left hand of Eqs (1) and (2), corresponding 
to the Arrhenius and Savata methods respectively and 
shown in the Methods section, are plotted vs. the 
inverse of the temperature considering several of the 
most usual models in the literature. The f(a) and g(a) 
functions are listed in Table 1. The resulting activa- 
tion energies, as obtained from the slope of the plots, 
and the regression coefficients showing the quality of 
the fits are included for the Arrhenius and Savata 



Table 1 f(a) and g(a) kinetic functions corresponding to the most widely employed kinetic models 


Mechanism 


Symbol 


f(a) 


f(a) 


Phase boundary controlled reaction (contracting area) 


R2 


2(1-a) ,/2 


2[1-(1 -a) V2 ] 


Phase boundary controlled reaction (contracting volume) 


R3 


3(1-a) 2/3 


3[1-(1 -a) 2/1 ] 


First order kinetics or Random nucleation followed by an 
instantaneous growth of nuclei. (Avrami-Erofeev eqn. n =1) 


Fl 


(1 - a) 


-ln(1 - a) 


Random nucleation and growth of nuclei through different nucleation 
and nucleus growth models. (Avrami-Erofeev eqn #1.) 


An 


n(1 - a) [-ln(1 -a)] 1_1/ " 


[-ln(1 -a)] ] - yn 


Two-dimensional diffusion 


D2 


1/[-ln(1 - a)] 


(1 -a)ln(1 -a) + Q 


Three-dimensional diffusion (Jander equation) 


D3 


3(1 -a) 2 ' 3 
2[l-(l-a)" 3 ] 


[1 -d -a) v ¥ 


Three-dimensional diffusion (Ginstling-Brounshtein equation) 


D4 


2[(l-a)-" 3 -l] 


(1 -2a/3)-(1 -a) 2/3 


Random Scission L=2 [10] 


L2 


2(a 1/2 -a) 


-2ln(a l/2 -1) 


Random Scission L>2 [10] 


Ln 


No symbolic solution 


No symbolic solution 
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Table 2 Activation energies and regression coefficients 
obtained from fitting the data from the simulated curve 
in Figure 1 to some of the most common ideal models 
employed in the literature, according to the Arrhenius 
method 



Model 


Corr. Factor r 


Ea (kJ mol" 1 ) 


F1 


1.000 


150 


R2 


0.980 


122 


R3 


0.992 


131 


A1.5 


0.999 


96 


A2 


0.999 


69 


A3 


0.999 


43 


D2 


0.982 


257 


D3 


0.989 


266 


L2 


0.993 


100 



The activation energy is determined form the slope of the plots in Figure 2. 



methods in Tables 2 and 3, respectively. Additionally, 
Figures 2 and 3 shows a selection of the plots resul- 
ting from the fitting of the data to the different 
kinetic models, so that the quality of the fit is clearly 
illustrated. Thus, when the simulated curve is 
analyzed by the Arrhenius method, six out of nine 
models deliver excellent fits to the data, with regres- 
sion coefficients over 0.99. Therefore, it is not really 
possible to effectively discern the correct model with 
this procedure. Moreover, as it can be observed by 
the values in Table 2, the activation energy obtained 
from the analysis is highly dependent on the kinetic 
model assumed, with only the fit to Fl yielding the 
correct value. Consequently, without further evidence 
regarding the correct kinetic model, the activation en- 
ergy cannot be established. The results are even more 

Table 3 Activation energies and regression coefficients 
obtained from fitting the data from the simulated curve 
in Figure 1 to some of the most common ideal models 
employed in the literature, according to the Savata 
method 



Model 


Corr. Factor r 


Ea (kJ mol 1 ) 


F1 


1.000 


153 


R2 


0.997 


141 


R3 


0.999 


145 


A1.5 


1.000 


102 


A2 


1.000 


// 


A3 


1.000 


51 


D2 


0.995 


274 


D3 


0.999 


290 


L2 


0.999 


101 


The activation 


energy is determined form the 


slope of the plots in Figure 3. 



problematic when the Savata model is employed 
(Table 3 and Figure 3). Using such method the fit to 
all models tested are excellent and no clear candidate 
can be appropriately selected. 

Consequently, neither the model nor the activation 
energy can be deduced from the use of a model-fitting 
method of kinetic analysis and a single non-isothermal 
curve. Note that this conclusion is reached after analy- 
zing simulated, error-free data. When using experimen- 
tal data it is probable that even a higher percentage 
of the tested models will adequately fit such data. 
Thus, in order to unambiguously determine the kinetic 
parameters by a model-fitting procedure, a set of 
experimental curves, each recorded under different 
heating schedules, must be employed [3,6]. Then, a 
set of four curves simulated assuming the aforemen- 
tioned kinetic parameters and heating rates of 1, 2, 
10 and 20 K min' 1 were analyzed simultaneously 
using the Arrhenius method. The resulting plots are 
shown in Figure 4. As it can be clearly noticed, only 
when the tested model is the right one, Fl in this 
case, are the plots positioned along a straight line. 
Therefore, using a set of curves recorded under different 
heating schedules, it is possible to unambiguously 
determine the kinetic model and, consequently, the 
activation energy. 

As a final note, it should be considered that these 
methods are proposed for single step reactions, which 
can be described by a single kinetic triplet. When more 
than one process is taking place, each of them is 
expected to be defined by a different triplet. The 
thermogravimetric curves included in Marini's work 
explicitly show a complex, multistep process and, 
therefore, such model-fitting methods cannot be 
employed. It is then recommended to resort to iso- 
conversional methods or attempt the deconvolution 
of the contributing steps in order to study them 
independently [3,7-9]. 

The uncertain results provided by the Wyden- 
Widmann method in Marini's work can be likewise 
explained. It is still a model-fitting method since it 
assumes the process is driven by an «th order kinetic 
model. As it has been reported elsewhere, an «th 
order model is described by a mathematical function 
that cannot replicate the initial induction period typ- 
ical of the function describing a chain scission model 
[10]. Thus, being the wrong model to describe the 
reaction, it is understandable that the required lin- 
earity is not achieved along the entire temperature 
range, as described in the paper. Had the experimen- 
tal data been fitted to the right model, such depend- 
ence of the activation energy on the number of data 
points considered would most probably have not 
been found. 
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Figure 2 Plots obtained from fitting the data from the simulated curve (E a =150 kJ/mol, A=10 s , and a F1 kinetic model) in Figure 1 
to some of the most usual kinetic models by means of Arrhenius model-fitting procedure. 



Conclusions 

It has been shown that the activation energy cannot 
be reliably determined by applying model-fitting 
methods of kinetic analysis to data obtained under 
non isothermal experimental conditions. Thus, the 
use of a set of curves recorded under different 
heating schedules is necessary, as recently recom- 
mended by the ICTAC Kinetics Committee [3]. It is 
also important to consider the nature of the reac- 
tion under study since only one step processes can 



be analyzed by this methodology. More complex or 
multiple step reactions require the use of isocon- 
versional methods or the deconvolution of the indi- 
vidual steps. 

Methods 

The simulated curves were constructed using a 
Runge-Kutta 4 th order numerical integration method 
by means of the Mathcad software (Mathsoft, Needham, 
MA, USA). 
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Figure 3 Plots obtained from fitting the data from the simulated curve (E a =150 kJ/mol, A=10 10 s"\ and a F1 kinetic model) in Figure 1 
to some of the most usual kinetic models by means of Savata model-fitting procedure. 



The Arrhenius method is based in the following 
equation: 

where da/dt is the reaction rate, f(a) is the kinetic 
model, A the preexponential factor, E the activation 
energy and T the temperature in Kelvin. The Savata 
method relies in the following equation, which is 



obtained after integrating the expression above and 
reordenating terms: 

E fAE\ 
log\g(°c)] = -0.4567— - 2.315 + %(^— J (2) 

where g(a) is the integral form of the kinetic model. 
For determining the activation energy, the left-hand 
side of Eqs (1) and (2) are plotted against the inverse 
of the temperature. The value of the activation energy 
is then deducted from the slope of such plot. 
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Figure 4 Plots obtained from fitting a set of simulated curves (E a =150 kJ/mol, A=10 1C s"\ and a F1 kinetic model) assuming heating 
rates of 1,2, 10 and 20 K min" 1 to some of the most usual kinetic models by means of Arrhenius model-fitting procedure. 
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